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Using a slave boson representation of multi-band Hubbard models, we in- 
vestigate a two-band model relevant to layered perovskites in the vicinity 
of half-filling. Beside the strong influence of the Hund's rule coupling, we 
obtain that the phase diagram separates into two regions: a weak to mod- 
erate coupling region where the effective mass is weakly renormalized, and 
a strong coupling regime where it is strongly renormalized. The transition 
between these two regimes is very sharp. It takes place in a (vanishingly) 
small domain. A ferromagnetic instability in only found in the strongly cor- 
related regime, and is triggered by the Hund's rule coupling. The results are 
compared to La-doped layered ruthenates. 

PACS numbers: 71.10.Fd, ll.SO.h, 74.70. Pq, 75.30.Kz. 



1. INTRODUCTION 

The richness of the phase diagram of the Ruddlesden-Popper series of 
the ruthenates makes it a fascinating topics of current solid state physics. 
Indeed single-layered ruthenates can be a superconductor (S^RuO^ as dis- 
covered by Maeno et alB or an anti-ferromagnetic insulator (Ca2Ru04), 
while multi-layered compounds may be ferromagnetic metals (SrRu03 and 
possibly Sr3Ru207) or an anti-ferromagnetic metal (CaRuOa). This rich- 
ness is particularly amazing when one observes that all these compounds 
are, from a theoretical point of view, remarkably similar. Indeed Ca and Sr 
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are divalent cations, the electron-electron interaction is strongest on the Ru 
atoms, which are common to all compounds, and the lattice structure does 
not drastically change when scanning through this family. Nevertheless these 
small changes must be responsible for the variety of the ground states. This 
thus calls for an accurate description of the electronic band structure, for in- 
stance on the level of the tight binding approximation. Note that Ca2RuO/i 
appears to be a Mott insulator below 300 K. Indeed, transport! and optical! 
measurements reveal a narrow gap of about 0.2 eV. It orders magnetically 
below Tjy =100 K, and therefore the insulating state may not be attributed 
to the magnetic ordering. 

Another striking feature is provided by La doping (or other trivalent 
cations) of both Sr- and Ca-based single layered compounds: in both cases 
a few percent doping turns these systems into ferromagnetsQ'Q. For the Ca- 
based compound, the Curie temperature is about 100 K, while substitutions 
of Ca with other trivalent dopants, such as Y, Ce or Pr result in ferromagnetic 
states with Curie temperatures varying from 120 to 200 K§'i. For the Sr- 
based compound, La substitution results into a ferromagnetic state, with a 
Curie temperature of 40 K at a concentration of 4 percent of dopantH. 

Instead of investigating the nature and the cause of the superconduc- 
tivity in Sr2Ru04, which is still a controversial matter@l§i0llll, we rather 
concentrate on the Fermi liquid parameters and the magnetic structure of 
the entire series. The paper is organized as follows: we first introduce a 
model meant to capture the main feature of this series on a minimal model 
level. It includes the details of the band structure, inspired by a proposition 
from Noce and Cuocdp which is compatible with the experimentally deter- 
mined Fermi surfaced and spectroscopic datalia. With minimal changes it 
can be applied to other perovskites such as vanadates and chromites. Sec- 
ondly we explain how we are dealing with the strong Coulomb and Hund's 
rule couplings. Thirdly we present numerical results pointing towards a pic- 
ture where both Ca-, and Sr-based ruthenates derive from a Mott insulator. 
The former system appears to be a Mott insulator at commensurate filling, 
while the second one should rather be seen as a self-doped Mott insulator. 
This is due to a partial filling of all sub-bands. 



2. THE MODEL 
2.1. Band Structure 



The band structure of Sr2Ru04 has been determined in the LDA by 
OguchS and SinghllL Later on it turned out that it can accurately be 
reproduced by a simple tight binding model by Noce and Cuocdla. The 
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bands which are crossing the Fermi energy are involving the three 4d t2o 
orbitals of Ru, which are hybridizing with the 2p orbitals of O. They are filled 
with 16 electrons per Ru02 unit. Since the RuC>6 octahedra are elongated 
along the z-direction, the d xy orbital energy level is lying lower than the other 
two. From the LDA calculation, the corresponding band is nearly filled. In 
the following it will be neglected, except for possibly contributing to self- 
doping effects. The & x (y)z orbitals hybridize only with the p x (j/) orbital of 
the neighboring O along the z-axis, and with the p z orbital along the x 
(y)-axis. On top, the latter two orbitals weakly couple to one another as 
well. It numerically turns out that the dispersion in the z-direction is small, 
and will therefore be neglected. The resulting dispersions are given by the 
eigenvalues of: 
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\a = (4*»4*>Plz>pL) fc<y » and a k = -4t 3 COS (%^) COS (^) , e k = 
2t 5 sin(^p), and f k = 2t 5 sin(^). Here -?f < k x (k y ) < ?f and a 
represents the lattice constant (the Ru-Ru distance). It is from now on set 
to unity. The LDA dispersions is best reproduced with the parameter set 
t 3 = 0.1 eV, t 5 = 0.85 eV, e p = -2.4 eV and e d = -0.9 eV. Integrating 
out the oxygen bands, and neglecting their frequency dependence, yields the 
effective model: 
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up to a constant energy shift. Here we introduced: a k = —At'smk x sm.k yi 
e~k = t cos k x , and fk = t cos k y . The effective hoppings t and t' are related 

to the original parameters by t = 2-^ and t = ta ( — I ■ Note that the 
two-dimensional character of the dispersions 
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with v = ±1, purely follows from the off-diagonal terms. These two eigen- 
values yield the dispersion one obtains in LDA calculationt3. The resulting 
density of states is symmetrical with respect to half band filling, where it 
has a minimum. The band width W is 1.2 eV. 
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2.2. Interaction Terms 

We have modeled the Coulomb interaction using a Hubbard U -term and 
Hund's rule coupling term J#-terms, yielding the local Hamiltonian: 

Hint — U ^ ^H,p,u^i,p' ,a 
i,o,p' <p 

+ Ui ""AtVj + ^E^I^Al. ( 4 ) 

i,p'^p i,p 

where a (p) is a spin (band) index. The relation between the coefficients 

U n = U + nJ H (5) 

holds for perfect cubic symmetp0'El, otherwise one should work with the 
matrix elements Fq, F%, and F$B. In the following we nevertheless stick to 
the minimal model using the relation Eq. (||). This amounts to deal with 
the parameters U, Jh, t, and t' (with t setting the energy scale). Since 
the ruthenates (Sri-zCa^RuC^ may be Mott insulators (depending on the 
value of x), perturbative calculations may be problematic. Here we resort 
to a slave boson approach@@0, which is the generalization of the Kotliar 
and Ruckenstein representation of the Hubbard model£3, to the multi-band 
model. At this stage one may wish to introduce a rotation invariant formu- 
lation in the spirit of Li et al. or Fresard and Wolfle's representationsE3'cil. 
In the one-band model, such a representation turned out to be necessary to 
obtain the correct contribution of the spin fluctuations to the specific heat 
on the one-loop WeltjjJ. It also appeared very useful when investigate incom- 
mensurate statesc3'c3c3. In the systems under consideration, perfect cubic 
symmetry is absent, and therefore rotational symmetry in spin space is not 
a requirement. Thus the local model Eq. (|j) can be viewed as a model in 
its own right, despite of lacking rotational invariance. Note that the tightly 
related Gutzwiller Approximationhas been suitably generalized to deal with 
rotationally symmetrical modelsta. Manifestly anti-ferromagnetic solutions 
are playing a role as well, as discussed by HasegawaS. This is however a less 
puzzling issue than the occurrence of ferromagnetism, and is therefore not 
investigated here. Since there is no van Hove singularity at half-filling, we 
expect that an anti-ferromagnetic instability can only occur above a finite 
(and substantial) U, as it happens in the t — t' — U modeS or in the one 
band model on the hexagonal latticed. Recently Hasegawalll showed that 
orbitally ordered states in this model seem to rather appear for negative Jh ■ 
In the slave boson approach to the multi-band model one introduces 
one slave boson field for each of the sixteen atomic configurations, and four 
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Fig. 1. Effective mass renormalization at p = 2 for Jh/U= (thick full 
line), 0.01. (thick dotted line), 0.02 (thick dashed line), 0.05 (thick dashed- 
dotted line), 0.1 (thin full line), 0.2. (thin dotted line), 0.3 (thin dashed 
line), 0.4 (thin dashed-dotted line) and 0.5 (thin dashed-dotted-dotted line). 
The circles indicate the location of the first order transition. 

auxiliary fermionic fields f a (a being a composite spin and band index) as : 

l«> = P ] Jl\ vac ) i 



la, a') = S ,flf,\vac) 

i / aa a 1 



(6) 
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On top of those, there is a boson e related to empty sites, and a boson q 
related to four-fold occupied sites. The original fermionic operators c\ a are 
expressed in terms of the auxiliary fields as c] 



zjafia witri z ^ given by Eq. 



(4-5) of Ref. 19, while the constraints linking the auxiliary fields are given 



by Eq. (6) of Ref. 19. Altogether, after having integrated out the fermions, 
we obtain the grand potential in saddle-point approximation as : 
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Fig. 2. Effective mass renormalization off half-filling for Jh/U= 0.5 and 
p = 2.005 (full line), 2.01. (dotted line), 2.02 (dashed line), 2.03 (dashed- 
dotted line), 2.04 (dashed-dotted-dotted line), and 2.05 (thin full line). Inset: 
Blow up of the "metallic" solutions with the same parameters. 
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In a paramagnetic or a ferromagnetic phase, the dispersions are given by : 
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where the dependence on a is only effective in the ferromagnetic phase. 
Though the resulting mean-field equations only correspond to the saddle- 
point of an action, the latter gets variationally controlled in the limit of large 
dimensions^. It would be highly desirable to reach the degree of accuracy 



Ferromagnetism in a Realistic Two-Band Model. 



2.05 
2.04 
2.03 

Q_ 

2.02 
2.01 



m f • 



□ * o 



,0 j: O 

.*....()■■•□" 

— ' L 

0.5 



0.7 



U/W 




Fig. 3. Coexistence regions of the paramagnetic phases for Jh/U= 0.5 (thin 
lines) and 0.3, (thick lines). The dashed lines locate the first order transition 
lines, the dotted lines U&, and the full lines U C 2- All lines are guide for the 
eyes, only. 



and reliability of the conservin g a pproximations developed by the Karlsruhe 
group for impurity modela23E3E§E3, but to date it appears difficult to extend 
them to lattice models. Nevertheless mean-field calculations turned out to 
compare well with numerical simulations!^, even when calculating the charge 
structure factorS 



3. NUMERICAL RESULTS 

The saddle-point equations have been solved on a 800 x 800 lattice, 
at a temperature T = t/1000. Since the electronic density in the ruthen- 
ates under study is p ~ 2, we neglected four- fold occupancies and empty 
configurations. This approximation is justified in the vicinity of the Mott 
transition. However it breaks down for densities above three (below one), 
and for weak coupling, where our results should be taken with care. In 
particular it results into having a finite effective mass renormalization in 
the non- interacting limit, which should not be the case. As obtained by 
Biinemann et al@tl, and Klejnberg and Spalek^l, the Hund's rule coupling 
has a strong influence on the Mott transition. While the latter is second 
order for Jh = and p = 2, or for any Jh for p = 1 or 3, it becomes first 
order for finite Jh at half-filling as shown on Fig. |l[ The metallic solution 
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Fig. 4. Effective mass renormalization at p = 2.1 for Jh/U= (full line), 
0.1 (dotted line), 0.2 (dashed line), 0.3 (dashed-dotted line), 0.4 (dashed- 
dotted-dotted line), and 0.5 (thin full line). 



of the saddle-point equations ceases to exist at a critical value U c %, without 
corresponding to a diverging effective mass. Moreover the effective mass 
is at most renormalized by a factor five for Jh/U > 0.01, in contrast to 
the one-band case. On top of this metallic solution, there is an insulating 
paramagnetic solution characterized by a vanishing value of all bosons, but 
d X zi,yzi and d X z\,yz\i anc ^ therefore a diverging effective mass (for finite Jh)- 
It extends down to U c \ = 0. We also note that U C 2 is slightly larger than 
U c , where the energy of the metallic and insulating solutions coincide. As a 
result the effective mass renormalization is even more modest in the metallic 
phase. 

Once the system is doped the situation changes little by little. For small 
electron doping, the first order transition remains but gradually vanishes 
with increasing electron concentration as shown on Fig. ^. The metallic 
solution is only modestly affected, except for that it allows for decreasing 
values of z in the vicinity of the Mott transition. The insulating solution 
becomes metallic under electron doping, and the truly insulating state is 
only found for integer fillings. However the effective mass renormalization 
remains very large, and accordingly the quasi-particle residue is small. Under 
these circumstances, the system may well be strongly influenced by other 
interaction terms, as reviewed by Vollhardt et al.0, or disorder effects. This 
may explain why many transition metal oxides remain insulating even upon 
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Fig. 5. Inverse magnetic susceptibility for Jh/U = 0.5 and p = 2 (full line), 
2.1 (dotted line), 2.2 (dashed line), and 2.3 (dashed-dotted line). 

substantial doping, such as Lai-xCa^VOaiil. 

When coming from large U, the "insulating" solution is found to exist 
down to a critical value U c i, which is smaller than U C 2- Beside the fact that 
U c i = at half-filling for any Jh, it is a complex function of density and Jh 
otherwise. It is shown on Fig. [| for Jh/U = 0.3 and 0.5; together with U c 
and U C 2- The coexistence region is seen to be largest for Jh/U = 0.3. It gets 
somewhat smaller for larger ratio of Jh/U, but its size drops abruptly for 
Jh/U — > 0. For larger doping the effective mass renormalization is getting a 
smooth function of U, as displayed on Fig. |||. Again it strongly depends on 
Jh, either in the vicinity of the Mott transition since its location depends 
on Jh, or when it is small, even for large U. Even though the curves are 
now smooth, one observes that the jumps occurring in the small doping 
regime are merely replaced by a very sharp drop which is, roughly speaking, 
separating a good metallic region from a nearly insulating one. 

We now switch to the magnetic susceptibility. The latter is obtained 
by evaluating the Gibbs free energy G(m) at mi = 0.1 and rri2 = 0.01, and 
expressing the magnetic susceptibility \ as : 



y = •^1^2 (9) 

X 2(G(mi) - G(m 2 )) K) 

As shown on Fig. [5] for Jh/U = 0.5 the spin susceptibility does not 
diverge at half-filling, regardless of the value of Jh, though it is strongly 
renormalized at the Mott transition, where the curve stops. In the "metallic" 



Raymond Presard and Mathieu Lamboley 



0.3 I 1 1 1 1 1 1 1 1 1 1 r 




Fig. 6. Inverse magnetic susceptibility for p = 2.3 and Jh/U = (full line), 
0.1 (dotted line), 0.2 (dashed line), 0.3 (dashed-dotted line), 0.4 (dashed- 
dotted-dotted line), and 0.5 (thin full line). 



regime, as defined above, the renormalization is slightly but increasingly 
suppressed with increasing electron doping. Nevertheless, in this metallic 
regime, there is no divergence of x> but merely a renormalization of it by 
at most a factor five. This renormalization is comparable to the one of the 
effective mass, though somewhat larger. 

As one might anticipate, the Hund's rule coupling has a strong influ- 
ence on the magnetic susceptibility, as shown on Fig. ^, where we display 
the inverse magnetic susceptibility as a function of the interaction strength 
for several values of Jh/U, at fixed density p = 2.3. For Jjj = 0, x 18 
barely affected by the interaction, and is not found to diverge for U up to 
10 W. Diverging susceptibilities start to appear once the Hund's rule cou- 
pling is finite. The critical value of U decreases with increasing ratio Jh /U, 
but always stays above U c . We therefore have no evidence for a weakly 
renormalized ferromagnetic system. This is in agreement with the one-band 
Hubbard model, where a ferromagnetic ground state onlv takes place for 
very large interaction strengths and moderate hole dopingo. 

The instabilities of the paramagnetic phases are collected on Fig. fj], for 
several values of Jh/U. The range of stability of the paramagnetic phase 
is seen to depend weakly on density for large ratio of Jh/U. In contrast, 
it may extend to large interaction strengths for Jh/U = 0.1. On top there 
is a strong asymmetry around p = 2.5, which is more a consequence of the 
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Fig. 7. Instability line towards ferromagnetism in the density-U plane for 
J H /U= 0.1 (full line), 0.2 (dotted line), 0.3 (dashed line), 0.4 (dashed-dotted 
line), and 0.5 (dashed-dotted-dotted line). The circles are indicating the 
corresponding C7 c 's at half-filling. Inset: Instability line (full line), and first 
order transition line (dotted line), for Jjj /U=0.5. 



difference between U c for p = 2 and p = 3 (see Ref . |I^) , than a consequence 
of neglecting four- fold occupancy. The choice m\ = 0.1 in Eq. @ results into 
a somewhat overestimated stability range for the paramagnetic phase, but 
numerical instabilities appear when working with smaller values for mi. As 
displayed in the inset of Fig. 0, the instability lines connect to the first order 
transition line separating two paramagnetic solutions, where the latter ends, 
within numerical accuracy. No ferromagnetic solution with magnetization 
mi has been found for very small doping and U > U c . 

When comparing this phase diagram to La-doped Ca2RuC>4, we see that 
a small amount of electron doping turns a Mott insulator into a ferromagnet, 
in agreement with experimenti. However the shape of the instability lines 
does not appear to be compatible with the experimental situation met in 
La-doped Sr2RuO,4. Indeed no instability line extends to the very vicinity 
of half-filling in the metallic regime. This point may be easily overcome by 
noticing that the band involving the d xy orbital is not completely filledH3E3, 
but only carries about 1.8 electrons per Ru atom. One should therefore 
concentrate on the vicinity of p = 2.2. But then, increasing the electron 
density from p = 2.2 at fixed ratio U/W does not allow for going from a 
paramagnetic state to a ferromagnetic one. 
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Fig. 8. Double occupancy in the same orbital for Jh/U = 0.5 and p = 2 
(full line), 2.1 (dotted line), 2.2 (dashed line), and 2.3. The arrows indicate 
the location of the corresponding ferromagnetic instabilities. 



In order to reconcile the theoretical and experimental data, one would 
therefore need to assume that La doping not only changes the electronic 
density, but also strongly influences (reduces) the band width. An argument 
pointing towards this direction is provided by a recent structural study by 



the ruthenates may easily vary to a large extend, even under a temperature 
change. There is little doubt that this holds for the energy levels and 
e p and the overlap integrals, and therefore for the band width. Recalling 
that the ionic radius of La 3+ and Ca 2+ are very similar, one may expect 
that La doping indeed induces a reduction of the band width, since Ca 
substitution does it, leading to a ferromagnetic instability. Another reason 
for the discrepancy may follow from considering the d xy band for self-doping 
only, instead of including it in its own right into the^ialculation from the 
outset. According to a recent work by Anisimov et alt^b, it indeed seems to 
play an important role. 

Since the four-fold occupancy turned out to be negligible, one may won- 
der whether further simplifications could be introduced without qualitatively 
altering the results. Indeed the determination of the saddle-point, even in 
the simplest paramagnetic approximation is not an easy task, and, except 
for the non- interacting and Jh = limits, quite careful numerical work is 
required. To address this question, we are displaying the less probable dou- 
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Fig. 9. Double occupancy in different orbitals with opposite spin projections 
for J H /U = 0.5 and p = 2 (full line), 2.1 (dotted line), 2.2 (dashed line), 
and 2.3. The arrows indicate the location of the ferromagnetic instabilities. 



ble occupancies on Fig. ^ and |j| After introducing : A 2 = d 2 z ^ xz , + d 2 z ^ , 
and G 2 = d%. z ^ . + d 2 z \ yz p we see that they behave quite differently as 
a function of U and Jh- Indeed the double occupancy in the same orbital 
(A 2 ), which is corresponding to the largest local interaction is suppressed 
monotonically under an increase of U. It nevertheless keeps a substantial 
value in the metallic regime, where it may therefore not be neglected. In 
contrast it becomes very small in the vicinity of the ferromagnetic instability, 
and presumably neglecting A is not going to have a strong influence on the 
determination of the instability line, especially in the large Jh regime. 

The double occupancy in different orbitals with S z = (G 2 ) behaves 
quite differently as shown on Fig. ||. An increase of U enhances it in about 
the entire metallic region. It then sharply drops in the vicinity of the Mott 
transition, but (except for the very vicinity of the Mott insulating state) its 
value remains quite sizeable, even for strong Hund's rule interaction and large 
U. We therefore conclude that even the instability line cannot be accurately 
obtained by setting G = 0, in contrast to A. That such an approximation, 
(which is corresponding to the limit Jh — * oo) will strongly influence the 
instability line is not particularly amazing, since the latter does depend on 
Jh- But this study shows that the dependence mostly arises through G, 
while the impact of A is far smaller. Therefore the precise value of the 
parameter Us in Eq. (||) is only marginally relevant. 
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4. CONCLUSION 

In summary we have investigated a two-band Hubbard model with a 
realistic band structure, where no van Hove singularity appear at, or in the 
vicinity of, half band filling. We have shown that the first order metal to 
insulator transition that occurs at half-filling first extends to small but finite 
doping in a strict sense, and second only weakly smears out for larger doping. 
This divides the phase diagram into a "metallic" region, where the effective 
mass is only weakly renormalized and no ferromagnetic instability is found, 
and an "insulating" region, where the effective mass is strongly renormalized, 
and ferromagnetism is found. On top, there is a surprisingly small region 
in-between, which is interpolating between these two behaviors. 
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